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ABSTRACT 

The occurrence of low-amplitude flux variations in blazars on hourly timescales, commonly 
known as microvariability, is still a widely debated subject in high-energy astrophysics. Several 
competing scenarios have been proposed to explain such occurrences, including various jet plasma 
instabilities leading to the formation of shocks, magnetic reconnection sites, and turbulence. In 
this letter we present the results of our detailed investigation of a prominent, five-hour-long 
optical microflare detected during recent WEBT campaign in 2014, March 2-6 targeting the 
blazar 0716-1-714. After separating the flaring component from the underlying base emission 
continuum of the blazar, we find that the microflare is highly polarized, with the polarization 
degree ~ (40 —60)%±(2—10)%, and the electric vector position angle ~ (10 — 20) deg±(l —8)deg 
slightly misaligned with respect to the position angle of the radio jet. The microflare evolution in 
the (Q, U) Stokes parameter space exhibits a looping behavior with a counter-clockwise rotation, 
meaning polarization degree decreasing with the flux (but higher in the flux decaying phase), and 
approximately stable polarization angle. The overall very high polarization degree of the flare, 
its symmetric flux rise and decay profiles, and also its structured evolution in the Q — U plane, 
all imply that the observed flux variation corresponds to a single emission region characterized 
by a highly ordered magnetic field. As discussed in the paper, a small-scale but strong shock 
propagating within the outflow, and compressing a disordered magnetic field component, provides 
a natural, though not unique, interpretation of our findings. 

Subject headings: acceleration of particles — polarization — radiation mechanisms: non-thermal — 
galaxies: active — BL Lacertae objects: individual (S5 0716+714) — galaxies: jets 
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1. Introduction 


Blazars are known for their intense non-thermal 
emission and pronounced variability across the 
electromagnetic spectrum, resulting from the ef¬ 
ficient energy dissipation taking place in the in¬ 
nermost regions of relativistic, magnetized, and 
non-stationary outflows — “jets” — produced 
by active supermassive black holes in the cen¬ 


ters of the evolved galaxies (e.g., Begelman et al. 


1984 Meier 2012). However, our understanding 
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of the exact physical conditions of the emitting 
plasma in blazar jets remains limited. The occur¬ 
rence of rapid and low-amplitude flux variations on 
hourly timescales, commonly known as microvari¬ 


ability, or intra-day/night variability (e.g., Wagner 
& Witzel [T995[ and references therein), provides 
additional challenge in this context, as the amount 
of relativistic beaming required by the plausible 
explanation for such phenomena is, in many cases, 
rather too extreme to be reconciled with the cur¬ 
rently favored models for the jet formation in ac¬ 
tive galactic nuclei (AGN). 

Several competing scenarios have been pro¬ 
posed to explain rapid, either low- or high- 
amplitude variability in blazar sources. Some 
include extrinsic causes, such as gravitational mi- 
crolensing ( Watson et al.||199'9 Webb et al.|2000 |; 
others involve intrinsie origin, including purely 
geometrical effects (e.g., the ‘light house effect’, 
Camenzind & Krockenberger 1992), or various 


plasma instabilities leading to the formation of 
shocks, magnetic reconnection sites, and turbu- 


lence (see the recent discussion in, e.g., 

|Narayan 

& Piran||2012 Subramanian et al.||2012| 

Marscher 

2014 Sironi et al. 2015 Saito et al. 2015). Since 


the radio-to-optical emission continuum of blazars 
(and BL Lacertae objects in particular) is known 
to be due to the synchrotron radiation by ultra- 
relativistic jet electrons, the temporal behavior of 
the optical polarization can be used as a power¬ 
ful tool for diagnosing the structure of the blazar 
emission zone and the source of their variability. 

For example, a significant polarization degree 
of the synchrotron flux indicates an anisotropic 
distribution of the jet magnetic field, which may 
be related to either a large-scale uniformity of the 


magnetic field lines (e.g., Begelman et al. 1984 
Lyutikov et al. 2005|), or a tangled, chaotic mag¬ 


netic field compressed or sheared by the flow (e.g., 
Laingl 1980[ 2002 Hughes et al.||198'5{ [Kollgaard 
et al. 1990[ Cawthorne fc Cobb||1990[ [Wardle ^ 
al. 1994 Nalewajko 2009). Moreover, the high 


duty cycle of the blazar polarization microvariabil¬ 
ity revealed by multi-frequency optical monitoring 
(> 50%, e.g., Andruchow et al.||2005 Villforth et 
al. 2009) indicates an origin linked to changes in 


the physical conditions of the jet plasma, rather 
than some external (to the jet) or purely geomet¬ 
rical effects. 

In this letter we present the results of our de- 
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tailed investigation of a particularly prominent 
and well-resolved optical microflare detected dur¬ 
ing the Whole Earth Blazar Telescope (WEBdJ^ 
campaign targeting S5 0716-1-714 in 2014, March 
2-6. During the campaign, high-quality, simul¬ 
taneous multi-band optical flux and polarization 
measurements of the source have been gathered. 
S5 0716-1-714 is one of the brightest blazars of 
the “intermediate-frequency-peaked BL Lac” type 
(IBL), exhibiting persistent activity in all wave¬ 
bands, including radio, optical, and y-rays (e.g.. 


Bhatta et al. 2013 Rani et al. 20151. In particu¬ 


lar, the source is known to show rapid optical fluc¬ 
tuations on the timescales of minutes and hours 


(Fan et al. 

2011 

Bhatta et al. 

2013 Wu et al. 

2014 

1; high and variable optica 

polarization de- 


gree of about > 30%, along with fast rotations of 
the position angle of the electric vector, have been 
also observed (e.g., Larionov et al. 2013). The 
microvariability duty cycle of S5 0716-1-714 is very 
high (~ 80% - 90%; |Webb|[2007| |Hu et al.|[M4t 
when compared with that of other blazars or other 
types of AGN (see Goyal et al. 20131. The selected 
source is therefore an ideal target to study polar¬ 
ization and spectra properties of the blazar optical 
microvariability. Below we discuss the detection of 
the very high polarization degree of the microflare 
component separated from the underlying, slowly- 
varying emission continuum in S5 0716-1-714, indi¬ 
cating a highly ordered magnetic field in the sub- 
region of the jet responsible for the production of 
the observed flux enhancement. 


2. The 2014 WEBT Campaign 

The simultaneous photo-polarimetric observa¬ 
tions of S5 0716-1-714 analyzed here were obtained 
as a part of the WEBT campaign which lasted 
for five days in 2014, March 2-6. At that time 
the source was in a low activity-level state, in¬ 
creasing slowly its optical flux after the histori¬ 
cal minimum from the end of 2013. The gath¬ 
ered dataset consists of high quality, high cadence 
multi-channel (BVRI) flux measurements from 
several observatories; on the two occasions, lasting 
for about 20 h (56719.44-56720.27 MJD; hereafter 
“Epoch 1”) and 18 h (56721.70-56722.44 MJD; 
“Epoch 2”), the high quality optical polarimet- 
ric observations were obtained from St.Petersburg 


(LX-200), Crimea (AZT-8), Elagstaff, AZ (1.8m 
Perkins) and Kanata (1.5 m) telescopes. Detailed 
analysis of the entire data set collected during 
the 2014 WEBT campaign, along with the dis¬ 
cussion on data acquisition and reduction, will be 
presented elsewhere (Bhatta et al. in prep.). 


3. Analysis and Modeling Results 


Eigure 1 presents the photo-polarimetric dataset 
gathered in the R-band for the Epoch 2 of the 2014 
WEBT campaign targeting S50716-I-714; the flux 
measurements in the remaining filters are more 
sparse, especially in the B-band during the epoch 
considered. As shown, a very prominent and well- 
resolved microflare with the approximately sym¬ 
metric and almost exponential flux rise and decay 
profiles — and a sharp peak in between — was ob¬ 
served in the period 79-85 h (from the campaign 
starting time at Oh; see the dashed vertical lines 
in the figure). The total observed intensity of the 
source varied at that time by ^ 10% in Tyar — 2 h. 
In order to extract the main characteristics of the 
flaring component, first we define the polarization 
degree PD, the position angle of the electric vector 
X, and the polarized flux PF through the Stokes 


parameters Q and U (see Rybicki & Lightman 


1986), 


PD 

X 

PF 


\/Q2 + [72 


F 



VQ^ + , ( 1 ) 


where F is the total flux of the source. 

Next we assume that the analyzed microflare 
constitutes a separate emission component super¬ 
imposed on the underlying, slowly variable “back¬ 
ground” provided by the base emission continuum 
of the blazar. Due to the linearly additive proper¬ 
ties of total flux and the Stokes Q and U intensi¬ 
ties, one therefore has 


F = Fq + Fi, 

Q = Qo F Qi : and 
U = Uo + Ui, (2) 

where the microflare and the base emission compo¬ 
nents have been denoted by indices “1” and “0”, 


thttp://www.to.astro.it/blazars/webt/ 
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Fig. 1.— Photo-polarimetric optical data col¬ 
lected for the blazar S5 0716-1-71 during the Epoch 
2 of the 2014 WEBT campaign. Subsequent pan¬ 
els from top to bottom present the total observed 
R-band flux F, the corresponding total polariza¬ 
tion degree PD, the position angle of the electric 
vector X, and the Stokes parameters Q and U. The 
dashed lines in the panels (a) and (d-e) denote the 
modeled “background” emission continuum. The 
period analyzed in this paper is marked by the two 
vertical dashed lines (between 79 h and 85 h after 
the start of the campaign). 


respectively. In the modeling of this particular 
microflare, all the background intensities {Fq, Qq, 
and ?7o) are found from fitting the data collected 
just before (79-80 h) and just after (84-85 h) the 
microflare, assuming that these may change slowly 
(linearly) with time (see dashed lines in the cor¬ 
responding panels of Eigure [^. Once the base 
intensities are found, we subtract them from the 
total intensities, to obtain Fi, Qi, and Ui, which 
further enable us to derive the basic parameters 
for the flaring component PDi, xi, and PFi, us¬ 
ing the standard relations given in equations 0 . 

In the above analysis, we use the data exclu¬ 


sively in the R filter, for which continuous flux 
measurements are available in both the intensity 
and polarization for the entire duration of Epoch 
2. To derive the errors for the base intensities, we 
take the square root of the average of the variances 
of the data points used for the background com¬ 
ponent fitting (i.e., just before and after the mi¬ 
croflare). This is a reasonable assumption indeed, 
since in the absence of a flaring component (which 
increases the signal to noise ratio thereby reducing 
the measurement uncertainty; see Howell 1989| ), 
the errors of the flux measurement should corre¬ 
spond to those of the slowly-varying background 
continuum. The errors in the derived quantities 
{PD, X, and PF), for both the base and the mi¬ 
croflare emissions components, are derived using 


the standard error propagation formulae (Beving- 
ton &: Robinson||2003 ). 


Thus derived basic parameters of the flaring 
component in S5 0716-1-71 are presented in Eig¬ 
ure As shown, the microflare is highly polar¬ 
ized, with PDi ^ (40 — 60)%±(2 — 10)%, espe¬ 
cially when compared with the slowly varying base 
emission for which PDq ~ 8% ± 0.01%. Interest¬ 
ingly, the electric vector position angle of the flar¬ 
ing component, xi ~ (10 ~ 20)deg±(l — 8) deg, 
is only slightly different from that characteriz¬ 
ing the underlying background component, xo ~ 
30 deg ±0.8 deg. It is important to note that, ac¬ 
cording to the high resolution radio image ob¬ 
tained on 2014 Eebruary 24 within the VLBA-BU- 
BLAZAIl[^ project, xo is close to the innermost 
(within 0.12 mas from the core) position angle of 
the jet (^45°), while xi corresponds to the posi¬ 
tion angle of the jet farther down from the core 
(~20°). 


In Figure we also present the evolution of 
the flaring event in the Stokes parameters plane 
Qi vs. Ui (upper panel). As shown, there is an 
indication for a looping behavior with a counter¬ 
clockwise rotation, i.e. with higher polarization 
degree occurring in the decaying phase of the mi¬ 
croflare. In general, the Q — U diagrams illus¬ 
trate the evolution of the polarized flux together 
with that of the polarization angle. In the par¬ 
ticular case presented here, the observed behavior 
implies therefore a consistent rate of PFi changes 
(given by the distance between the consecutive 


thttps : //www .bu. edu/blazars/VLBA_GLAST/0716 .html/ 


4 



























0.4 


0.3 


0.2 


-- 6 -- ' 


0 0.1 0.2 0.3 0.4 0.5 0.6 

Q, (mjy) 


90 


70 


50 


•• • 


79 80 81 82 83 84 85 

Time (hr) 


t,, 

'""i. 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

F, (mJy) 


Fig. 2.— The derived basic parameters of the flar¬ 
ing component in S5 0716-1-71. Subsequent panels 
from top to bottom present the R-band flux 
the corresponding polarized flux PFi, the total 
polarization degree PDi, and the position angle 
of the electric vector xi ■ Dashes in the panels (c) 
and (d) denote the modeled “background” emis¬ 
sion continuum. The duration of the microflare is 
marked by the two vertical dashed lines. 


data points), with no significant rotation of xi- 
This is a very interesting result as some of the 
similar studies on optical variability show “ran¬ 
dom walk”-type behavior in the Q — U plane; 
whereas only few cases of them reveal a struc¬ 
tured evolution (e.g., Uemura et al.|2010 and ref¬ 
erences therein). Moreover, during the analyzed 
microflare the polarization degree of the flaring 
component decreases with the flux: as shown in 
the lower panel of Figure even though the un¬ 
certainties in the polarized flux measurements are 
large in the begining and at the end when flar- 


Fig. 3.— The evolution of the analyzed microflare 
in S5 0716-1-71 in the Stokes parameters plane Qi 
vs. Ui {upper panel), and in the Fi — PDi plane 
{lower panel). 


ing component nearly equals the baseline emis¬ 
sion, there is a clear counter-clockwise looping in 
the PFi vs. Fi plane, with the polarization de¬ 
gree reaching a minimum at the level of ~ 40% at 
the peak of the flare. We will come back to this 
interesting finding in the discussion section below. 

Finally, we attempt a limited spectral analysis 
for the microflare, using the available BVRI data, 
and considering the multi-band flux measurements 
to be simultaneous if collected within 10 min win¬ 
dows. After separating the base components of the 
microflare in all the filters, we evaluate the spec¬ 
tral indices of the base component before and after 
the microflare, assuming a power-law distribution 
oc In this way we find that the spectral 

shape of the base emission continuum is nearly 
constant in time (for the period analyzed), with 
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tto ~ 1.44 ± 0.01. Next, we estimate the spectral 
index for the flaring component from the resid¬ 
ual (total minus base) multi-band fluxes. Rather 
sparse BVI data allow us, however, to estimate 
only roughly the spectral slopes of the flaring com¬ 
ponent in the initial and final phases of the flux 
enhancement. The resulting spectral indices varies 
between ai ~ 1.15 ±0.07 and ~ 1.46 ±0.16; large 
uncertainties precludes us however from making 
any definite statements on the spectral evolution 
of the flaring component. 


4. Discussion 


A few studies exist in the literature where at¬ 
tempts have been made to separate a base emis¬ 
sion component from a variable emission compo¬ 
nent in the total and polarized flux optical obser¬ 
vations of blazars. In most of the cases, the au¬ 
thors assumed constant background, and derived 
relatively high values for the polarization degree 
of flaring components 20% — 50%), with the 
electric vector position angles typically aligned to 


the jet directions ( 

Hagen-Thorn et al. 

2008 

Saki- 

moto et ai.||2013| Morozova et al. 

2014 

Covino 

et al. 2015). Particularly interesting for the work 


presented here are the results reported by [Sasada] 
et al. (20081 for S5 0716±714 based on the photo- 


polarimetric Kanata data gathered on 2007 Oc¬ 
tober 20; these authors allowed for slow (linear) 
changes in the base emission continuum, deriving 
for the variable component PD « 27%, with y « 
150 deg basically constant during the flaring event. 
For the flare analyzed in this paper, we found a 
much higher polarization degree of ~ 40% — 60%, 
and the polarization angle ~ 10deg—25deg only 
slightly misaligned (by < 20 deg, at most) with 
respect to the polarization angle of the base com¬ 
ponent, or to the position angle of the mas-scale 
radio jet. 


The observed high polarization degree of the 
flare, its symmetric flux rise and decay profiles, 
and also its structured evolution in the Q — U 
plane, all imply that the observed flux variation 
corresponds to a single and well-defined emission 
region characterized by a highly ordered magnetic 
field. A small-scale but strong shock wave prop¬ 
agating within the outflow, and compressing effi¬ 
ciently a disordered small-scale jet magnetic field 
component, is a natural and often invoked inter¬ 



e [deg] 


Fig. 4.— The expected polarization degree PD 
from a compression at the front of a strong shock 
as a function of the angle between the shock nor¬ 
mal and the line of sight 8, calculated using the 
equation ([^ for three values of the bulk Lorentz 
factor F = 10, 15, and 20 (dotted, dashed, and 
solid curves, respectively) with a = 1.5 and 1.0 
(thick and thin curves, respectively), and fc = 1/3 
and 1/2 (black and blue curves, respectively). 


pretation. And in fact, the stable electric polariza¬ 
tion angle positioned along the jet axis can easily 
be reconciled with the shock interpretation. More¬ 
over, in the case of a propagating shock, one may 
naturally expect a counter-clockwise looping be¬ 
havior in the F — PD plane with the polarization 
degree correlated with the flux, resulting from the 
increasing shock compression in the initial phase, 
followed by the decompression and spectral steep¬ 
ening due to radiative cooling of the emitting par¬ 
ticles at the later stages of the shock evolution 
(e.g., Hagen-Thorn et al.|2008 see also in this con¬ 
text the discussion in Perlman et al.|[^ll| regard¬ 
ing the optical polarimetric data for the flaring 
HST-1 knot in the M87 jet). In the case of the 
dataset analyzed here we do observe the counter¬ 
clockwise looping, however with the polarization 
degree anti-correlated with the flux. 


In order to discuss this finding in more detail, 
though still only qualitatively, let us first note 
that the polarization degree expected from a shock 
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compression can be expressed as 

PD = (1 - fc^) sin^g 

5 + 3q! 2 — <5^ (1 — fc^) sin^ 6 ’ 


where a is the spectral slope of the non-thermal 
emission continuum, k is the degree of a unit 
length shock compression, and <5 = l/r(l — 
P cos 9) is the Doppler factor of the emitting 
plasma characterized by the bulk Lorentz factor 
r = (1 — and the angle between the shock 


normal and the line of sight 9 (e.g., Kollgaard et 
aL]|1990 ). Here we evaluate polarization degree as 


a funtion of 9 for a number of parameters rela- 
vant to our study. Following our spectral analysis 
presented at the end of §[^ we chose the two lim¬ 
iting values of the spectral index a = 1.0 and 1.5, 
and set fc = 1/3 as appropriate for a strong and 
relativistic (in the jet rest frame) shock; we also 
consider a weaker shocks with k = 1/2 for com¬ 
parison. We then take the three representative 
values of r = 10, 15 and 20, with the limitation of 
9 < 10 deg, for rather illustrative purposes only, 
noting that the analysis and modeling of the radio 
interferometric data regarding structural changes 
in the mas-scale jet of S5 0716-1-71 implies the jet 
viewing angle 9j < 5 deg and the jet Doppler fac¬ 


tor 10 < 5j < 30 (e.g., Bach et al. 2005 Rani et 


^|M5 and references therein). 


With the given free parameters, we find that 
the expected value of the polarization degree de¬ 
pends strongly on the combination of F and 9. 
As illustrated in Figure even small changes in 
both parameters may result in significant changes 
in PD whereas the change in a by 0.5 may not 
produce appreciable change in PD. This exercise 
suggests in particular that the observed counter¬ 
clockwise looping behavior with the polarization 
degree anti-correlated with the flux could possibly 
be explained by assuming that at the initial stages 
of the shock evolution the shock normal starts to 
deviate from the jet axis (in a sense 9 > 9j), and 
that after the peak of the flare the disturbance 
decelerates, aligning finally its direction to that 
of the large-scale outflow. Interestingly, a care¬ 
ful look at the xi evolution in the bottom panel 
of Figure suggests that such a scenario may in¬ 
deed be the case. Also, the maximum polarization 
degree expected in the model for the maximum 
shock compression corresponding to fc = 1/3, is 
~ 60%, which is nicely consistent with the upper 


bound for PDi derived in §[^ Finally, we note 
that in the framework of the above interpretation 
one may expect the other microflares in the source 
to be characterized by substantially lower polar¬ 
ization degrees and larger misalignments between 
the polarization vector and the jet position an¬ 
gle (when compared with the microflare analyzed 
here), due to the inevitable spread in the kine¬ 
matic parameters F and 9 of small-scale shocks 
developing within the outflow. A similar geomet¬ 
rical approach for the interpretation of flux and 


PD variability was adopted by Raiteri et al. (2012) 


and Raiteri et al. (2013) while analysing data from 
long-term WEBT observations, to account for the 
correlation between flux and PD in the blazar 4C 
38.41 and anti-correlation in BL Lacertae. 

The above interpretation is obviously not a 
unique one, and several other models — including 
for example a helical distortion in strongly mag¬ 
netized plasma subjected to MHD instabilities — 
may possibly account for the observational find¬ 
ings as well. Also, if the flaring zone consists of an 
underlying, highly uniform longitudinal magnetic 
field on top of which a shock propagates injecting 
freshly accelerated particles, it is plausible that 
the observed F — PD anti-correlation is due to a 
reduction in the net polarization by an increasing 
transverse field related to the shock compression 
of a random, small-scale magnetic field component 
(see in this context Cawthorne et al.]|1993 ). Only 
a regular, multi-band photo-polarimetric moni¬ 
toring of the source with the time resolution of 
the sub-hour scale, allowing for a more precise 
characterization of multiple optical microflares in 
S5 0716-1-714, could tighten modeling constraints 
and disprove various alternative scenarios, provid¬ 
ing in this way a unique insight into the small-scale 
structure of relativistic outflows in AGN in gen¬ 
eral. What is the robust conclusion from the anal¬ 
ysis presented in this paper, however, is that the 
small-amplitude flux changes observed in blazar 
sources are related to uniform, coherent emission 
regions characterized by a highly ordered mag¬ 
netic field and small linear sizes of the order of 
I ~ 10^® (T^ar/hr) ((5/10) cm, and as such consti¬ 
tuting (most likely) only small sub-volumes of the 
outflows. 
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